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HfiCluMl and davice for woaBuriao* detaxminiog and controlling 
flatnciaa e£ a metal strip 

5 veolmieal riald o£ tha Invention 

The invention relates to a method and a device for measuring, 
determining and controlling flatness of a substantially long 
and flat sheet or strip of material, such as copper, steel, 
10 aluminium or other materials such as plastics, in a continuous 
or semi -continuous production process. 

More particularly it is a method and device for measuring, 
determining and controlling flatness in a cold or hot rolling 
15 mill . 

The invention further deals with a method and device for 
measuring, determining and controlling flatness using human- 
machine interfaces and control algorithms. 

20 

Baokgroand o£ the invention 

In the rolling of strip and sheet materials it is common 
2S practice to roll a material to desired dimensions in a cold or 
hot rolling mill stand and then feed the resulting strip to a 
coiler where it is wound up into a coil. In figure 1 -is- a 
flatness control is illustrated that typically comprises of 
three main parts, measurement, control and actuators (that 
30 correct the flatness error) . tPhe difficulty with the flatness 
control consists on achieving a very good flatness. The reason 
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why the flatness is not alws^ys *«• perfect, is that all three 
parts have a number of limitations. 

The measurement is for example often iSiscrete both across the 
5 width of the strip, as it has a limited number of measurement 
zones, and along the strip, as there are only a limited number 
of ^e«nF measurements done per rotation {say, four measuring 
zones in the rotational direction) . The control has a complex 
mathematics that cannot always he performed in real time and 

10 the actuators are only corecting certain type of errors and it 
is not sure they can correct any shape of flatness error. For 
example, a bending and a tilting actuator can correct only a 
shape that is a sum of a parable and a straight line, but of 
course the flatness error can in principle have any irregular 

15 shape. 

A metal strip that is subject: to different degrees of reductioj 
across its width will be elongated in varying length over diff- 
erent sections of the strip. Normally/ the reduction results ii 
20 certain pre-defined stress profiles for the strip. 

The difference between the measured strip flatness and a targel 
flatness curve is defined as the "flatness error* . Practically 
the flatness error of the strip can be gauged by a special 
25 measuring roll, e.g. a Stressometer roll, an ABB product. 

Prior Art 



30 



Several different methods have been developed to correct the 
flatness error by a so-called flatness control operation. All 
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of these methods are based upon local modifications of the gap 
between rolls. 

The tension in the strip between a mill stand and a coiler is 
5 carefully monitored and it is known to measure tension distri- 
bution across a strip in order to regulate the flatness of the 
rolled material. 

In US 3,481,194 Sivilotti and Carlsson disclose a strip flat- 
10 ness sensor. It comprises a measuring roll over which the strij 
passes between a mill stand and, for this example^ a coiler. 
The measuring roll detects the pressure in a strip at several 
points across the width of the strip. The pressure represents « 
measure of the tension in the strip. Ilie measurements of 
IS tension in the strip result in a map of flatness in each of 
several zones across the width of the strip, 

US 4 r 400, 957 discloses a strip or sheet mill in which tensile 
stress distribution is measured to characterise flatness* The 

20 measures of flatness are coxnpared to a target flatness and a 
difference between measured flatness and target flatness is 
calculated, as a flatness error. The flatness error is fed 
back to a control unit of the mill stand, so as to regulate an> 
control flatness in the strip in order to approach a zero 

25 flatness error. 

Discrete flatness sensors and rolling mills with several 
actuators have drawbacks in that the sensors are not measuring 
continuously* This means that certain types of flatness errors 
30 are not detected, in addition to this, due to physical 

limitations, actuators cannot correct any type of flatness 
error, even if these are detected. 
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One further disadvantage is increased complexity in sensor 
processing, actuator response prediction and the need of 
precise mathematical mappings between these sensors and 
5 actuators. 

Normally the optimization of the flatness is done by different 
coxt<pensations of the space inapping function. These involve a 
great deal of numerical and process control methods. Typical 
10 examples are artificial neural networks, multivariable control, 
adaptive control, etc. These are difficult for plant operators 
to understand and handle. 

Moreover* though additional information nay be available for 
15 the operator from e.g. visual inspection of the strip, the 
operator cannot easily translate this information into the 
existing mathematical formalism of the flatness control. For 
example, the operator may see certain patterns on the strip 
like "quarter-buckles", taut it is not obvious how he may use 
20 this information to amend the flatness control paradigm. 

•A 

In the following the background of control algorithms for 
; .-. control equipments will be described for the purpose of later 

• * M 

: understanding of the inv^ition. In doing this, roforona 

25 references will be made to the following parameters: 

: :'t N = Nuiriber of raeaeuranent zones (nittriber of samples for actuato. 

r : function) 

P = Number of actuacore 
I : 30 Jbi =: Flatness error for zone 1 (i=l,..MN) it is the flatness 

error space 
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ax^i ■ Actuator function at maasureraant zone i for actuator j , 

A = The mill matrix, with components at,j used for space mapping 
Xj s actuator parameter to calculate, normalized in range 
5 [0, . . ,1) f it is the actuator space 

s K&Boia the transpose of a matrix A 
R" = It is the space of real n-dimension vectors {Xn xa,-, Xn) . 

The term space denotes a vector space. We use the vector space 
10 o£ matrices o£ different sizes and the vector space of 
polynomials . 

The mill matrix A is obtained from the mill manufacturer, from 
measurements or identification aerations, and the flatness 
15 error matrix Jb is measured on-line with the measurement roll. 
These matrices are related to the actuator parameter matrix x 
by the following formula: 



20 



Ax-b (1) 



Flatness control means generally finding the vector x when the 
vector A is a given constant and b is measured on-line. 
Since A is not a square matrix (in general N > P) , there are 
more eqcaations N than tmknowns P for the system of eqpjiations 
25 represented hy (1) . That means, in general, there is no a 
unique solution but there can be found an x such that Jb is 
"close* to Ax. Here ♦'close" is in tewns of a Euclidian distanc 
expressed as a norm in the space. 

30 If X varies and if the columns of A are linearly independent, 
then the product Ax describes a space on called ''the colmm 
— » a* Y-anna nf A« - notation R(A) , or *the span 
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of the columns of A*. It is knovm from elementary geometry in 
2D/3D spaces and from functional analysis in N-space that the 
closest distance is given by an orthogonal projections the 
smallest distance between b and Ax is the projection of b on 
5 the apace of Ax, see figure 2. 

Two arbitrary vectors v and w are orthogonal when their inner 
product denoted <v, w> is aero. For a matrix vector space, the 
inner product is defined as the product vV^ so orthogonality 
10 means v'^w=0. Por our case, let r«i-Ajcbe the residual {the erro 
between a solution Ax and the vector b) > The smallest residual 
is the orthogonal projection of r onto Ax, which according to 
the definiticm is: 

15 Q^A^r^A^ib-Ax) (2) 

By separating the terms, we obtain: 

A^Ax^A^b (3) 



20 



This is called in literature the jaormal eguation. The factor fi^ 

is called the Gram matrix of A (is always a symmetric 
square matrix) . Multiplying {2) with (A^A)'^ we obtain x: 



. 25 x^(A^AyWb=A'^'-b (4) 



The term A*^ ^CA^Ay'A^ ±6 called the left pseudo-inverss matrix 
of A. 

30 The formula (4) above is used for the ^*model -based method*' in 
the existing Stressometer 5.0/6.0 flatness control etjuipment 
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and also by a majority of other producers of flatness control 
equipment. Other descriptions of this method appear under the 
names such as least square method" i orthogonal vectors" or 
**residual square minimization*, however, all of these are less 
5 general formulations of the vector space and orthogonal 

projection described above- For exairrple orthogonal vectors may 
be obtained by a Qram-Schmidt ortogonalissation procedure; 
however the ortogonalization is just another way of computing 
the pseudo- inverse. 

10 

Jmsgre morphinfir {short for image metamorphosis) is a popular 
technique for creating a smooth transition between two 2D 
images, Morph is the process of changing one object to another 
one. 

IS 

In Amorphium, for example ^ one can morph organic shape 1% to 
100%, Morphing organic shapes or shapes with the same nuniber oi 
polygons and vertices is called basic morph. All of Amorphium 
basic object (organic) have the same nuiriber of polygons and 
20 vertices and therefor can be morphed 100%. 

Morphing synthetic or shapes with different number polygons an< 
vertices is called shape morph. With shape morph Amorphium 
attempts to build a map between the vertices of one object and 

25 the vertices of the other- Therefor you do not end up with the 
morph shape looking identical to the one you morph it with. Al 
the standard shapes are synthetic ^ so are imported shapes 
FontMan, Keshman and shapes that have been merged. But if you 
morph an object that is synthetic, with an altered duplicate o 

30 it, it will morph 100% because it is considered basic morph. X 
has the same number of polygons and vertices. 
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The method can be applied to volume-based representations of 
objects and the result is a snvooth transition between two 3i3 
objects (3D metamorphosis or 3D inoxphiag) . There are several 
different algorithms to create the transition depending on the 
intended function. 



Sumnaxy of the Invention 

10 As Strip products with a high degree of flatness are essential 
for the conpetitiveness of mills today, the object of the 
present invention is to solve above indicated problems and 
present a new method and a device for measuring, determining 
and controlling flatness of a substantially long and flat sheet 

15 or strip of material. 

Flatness Control, PC, essentially maps points moving in the 
multidimensional space of the measured flatness error into 
points decemined in the multidimensional space of the actua- 
20 tors such chat certain optimality conditions are met. This is 
the traditional mapping used today in PC, we call this as the 
'first mapping' . 

The first mapping is often difficult to optimize due to the 
2S nature of the actuator space and due to the sensor data vs. 
reference values used for the eospputation of the flatness 
error. For instance the coir^lece action of all the actuators i 
not linear in the individual actuators, it cannot produce the 
optimal mapping or the mapping {possible conibinations of 
30 actuator positions) is not unique. 
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Another factor that limits the control is the discrete nature 
of the measurement. A discretization is inherent with most of 
the current: measurement technologies both across the strip 
width (the weasuremenc jsrones of the measuring roll) and along 
5 the strip length (sensors measuring when in contact with the 
strip) . 

The present invention proposes a solution to the problems 
discussed above by introducing a second mapping. This second 

10 mapping can be made from reference models such as video samples 
or mathematical functions to output space optimization 
matrices. The reference models enter into the computation of 
the main mapping. The result is increased accuracy and locality 
of the flatness control, A further advantage is a vast increase 

15 in sampling granularity that translates into more accurate and 
timely flatness control. 

The method of the invention is achieved byi 

- that a second mapping is made to optimization algorithms 
^20 acting in actuator model and sensor spaces using visual or 

image information, and 

- that the measured flatness vector is enhanced with 
information from a virtual strip. 



25 The present invention apply to any type of flatness measuremeni 
that is discrete. 

The present invention creates a mathematical map between the 
sensor space and the actuator space that improve substantially 
30 the deficiencies with known controlling methods. This comes 

from that extra information is entered into the system from th 
user^ or from a visual sensor, about the nature of the flatnee 
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error and that this information is used in a better way to 
optimize the actuator behavior, for any type of known 
controller. Acordingly, a visual equipment and visual 
operation, on sensor data, such as sensor fusion or inorpbing is 
used in the method according to the invention, as well as a 
control architecture that implements a space mapping. 

A method according to the present invention will increase the 
control speed substantially as the response time decreases. 



8ri«£ Daaeripticm oC nrasriagB 

For better understanding of the present invention, reference 
IS will be made to the below drawings /figures. 

Figure 1 illustrates a known flatness control equipment. 
Figure 2 illustrates a simple orthogonal projection. 
Figure 3 illustrates a control system according to the 
20 invention. 

Figure 4 illustrates a mathematical model generated strip. 
Figure 5 illustrates a number of examples of flatness profiles 
Figure 6 illustrates the fine-tuning of the reference model or 
the synchronisatifflci of the samplixig* 
25 Figure 7 illustrates an exanqple of a continuous strip 
prediction. 



30 



Detailed Daaeriptloa of Preferred Eabodlmenns 

Flatness control, as used in metal rolling mill industry is, 
from a mathematical perspective, a mapping of an abstract spac 
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of sensor-originated vectors into a space of actuators- When 
fine-tuning the control to specific mill conditions, one uses 
traditional paradigms corresponding to the iniplemented 
mathematical level of the control/ such as ''transfer function* r 
5 ^'poles and zeroes in the. state space", '^dynamic behavior'', 
"least-squares" and different parameterization or penalty 
functions. With the advent of the IT and virtual reality 
architectures^ it is increasingly the case that virtual strip 
models are available from video captures or mathematical 
10 models. Thus the idea is to move the tuning operation from the 
control abstract level up to the real physical properties that 
can be seen by an operator when using different visualisation 
and virtual reality techniques, 

15 The present invention presents a method to use the information 
contained in the virtual strip models to iiriprove the flatness 
control algorithm. The main features of the invention are: 

• The measured flatness vector is enhanced with information 
from the virtual strip using a morphing technique 
{enhancement of the first space mapping) . 

• By this operation the sampling granularity increases. That 
result in better and more accurate flatness control. 

• Since sampled data is available in shorter time interval, th 
control speed increases (response time decreases) 
substantially, 

• The virtual strip is mapped to a vector space optimization 
algorithm- This in effect is a second vector mapping that 
performs a fine^ local adjustment of the control optimizaeic 
(creation of a second space mapping) . 
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m short, the invention states that increased flatness control 
precision can be achieved by a two-way space nuapping in which 
visual or other continuous or near continuous a-priori strip 
information enhances <a) the measurement (b vector) and a 
5 traditional control loop and {b) the optimization algorithms 
that tells how x is computed such that to minimiae the rezidual 
lAx - b| . 

In figure 3 is illustrated a control system according to the 
10 invention, where: 

• A is a measurement roll, that produces a vector of flatness 
values per each measurement zone. 

• B is a Morphing unit (software in existing equipment or 
dedicated device) , that is a part of the HMI poftwaro . 

15 • C denotes reference strip models used for morphing- These ari 
picture files generated from video or using laser captures o; 
results from mathematical modelling. 

• D are target stress profiles (reference) . For the purpose of 
the method, the granularity of the target stress profile is 

20 the same as the granularity of the reference strip models. 

• E is flatness control computer uolng o.gt executing the 
algorithm described above. 

• P is are space conversion matrices that describe 
machematically the local properties of the strip (edge, 

25 middle etc) . Each matrix is mapped to a reference strip 
model. 

• G is a mill model, that is part of the existing flatness 
control, and it is a matrix («the Mill Matrix") . 

• H are actuators that correct flatoess. 



30 



The mill operator has excra information about the nature of th 
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inspection. By contrast, the measured value contains only a 
truncated part of all this information due to the discrete 
sainpling. A video frame can contain informtion at sub-milli- 
mster granularity while the flatness measurement may have only 
granularity between XO cm and up to 0.25 m. For example a video 
camera (or visual inspection) can show clearly so-called 
-quarter buckles" but when sampling with a measurement roll, 
these buckles are not completely detectable. This is 
illustrated in Figure 6: on a strip picture generated from a 
video capture, discrete roll measurement samples may fall 
anywdxere inside or outside buckles such that the buckles cannot 
be deduced from roll measurements alone. The visual infomatior 
can be generated also from mathematical models (so called grey, 
white or black models), an example of a mathematically 
generated model strip is illustrated in figure 4, According to 
the invention all types of generated or photo/video captured 
images can be used. 

According to the present invention, this information is added 
to the measured information from the reference model in order 
to improve the space mapping optimization. The information is 
typically roorphed. 

The operator has a number of reference models that have as 
output stress profiles, see figure that can be roorphed with 
the real measurement data. These models are obtained for 
exaitple from digital translations of video captures. 



To accomplish the morphing the selected reference model need 
be scaled and synchronized to match the measured data- 
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one alternative to do this is to make a visual mark on the 
strip (from video capture) where the measuretnent was done. This 
can be done for example by using a laser ray. 

5 in the graphical user interface the operator (or a pattern 
matching algorithm) can fine-tune the reference model hori- 
zontally to match the measured data. This method is shown in 
figure 7. 

10 The morphing algorithm modifies the reference model according 
to the measurement data resulting in a data set that is closer 
to the reality. 

The morphing function ean be calculated in different ways but 
15 essentially every measured data point affects several reference 
model poincs gradually. 

This new data set can be used to calculate a flatness error 
that is more accurate as compared to the existing system, 

20 The target profile is subtracted from the morphed profile (a 
vector subtraction) in order to compute a flatness error. This 
operation can be performed if the granularity of the target anc 
the morphed signal is the same. High granularity target 
profiles can be obtained by singly filling in r^etitively the 

25 same profile until the regaired granularity is achieved, or by 
using pictures or video captures from a de post-facto strip 
processing as described in the Patent EP 1 110 635 Al Method 
and device for controlling flatness (Jonsson, Meyer, 1999) . 
The resulting flatness error has a much higher granularity thai 

30 what is traditionally obtained from the original measured 

samples. This has the effect that the flatness control will be 
more accurate, faster and has a better prediction capacity. 
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The mill operator has extra information about the nature of the 
flatness fault, e.g. from a video camera or visual inspection. 
By contrast, the measured value does not contain all this 
5 inform tion due to satcQpling. For exanple a video camera {or 
visual inspection) can show so-called ^'quarter buckles* but 
when sanqpling with a measurement roll, these buckles are nob 
completely detectable. 

10 l*he operator has a number k of reference models of pre-defined 
stress profiles* These models are obtained from for exairgple 
digital translations of video captures. 

Information about the nature of the flatness fault can be used 
15 to design the flatness control algorithm to concentrate on 
coR^ensating zones where the problem is most significant, if 
for exanqple the strip has the pattern ^*wave edges' the 
actuators need to compensate strip edges rather than the 
center . 

20 

The present invention propose that there is established a one- 
to-one mapping: each reference model of the pre-defined stress < 
profiles has associated a matrix that describe the weight of 
each measurement zone for the flatness error correction. These 
25 weight matrices are used in the flatness control algorithm. 

The operator, or a pattern matching software, decides witch 
reference model is best to use. The reference model is mapped 
to one of the weight matrices. This mapping is the second 
30 vector space mappinff, and it is a key element of this 
invention. 
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There are some patents from e.g. Hitachi/Siemens, that use 
certain pattern matching techniques. However ^ these patents are 
focused on the use of certain algorithms {fuzzy, neuronal nets) 
for enhancing measureipent and are not refering to the control 
algorithm as this invention does. In other words, these patents 
claim improvements in what we call the fir^t space mapping. By 
contrast/ the present invention claims that substantially 
increased flatness control accuracy can be achieved by the dual 
space mapping. 



Mathematically, that means each reference model Ma (where i is 
the index of the model} « is znapped to a unique matrix Ki. 7he 
matrix Ki is a square matrix of size N x N (N being the nixmber 
of measurement zones) . The matrix has zero in all positions 
15 except the main diagonal, where there is the weight desired fo£ 
the action of actuators on the particular zone of size. 

The near-continuous space obtained after morphing^ gives, aftei 
a measurement sampling, a degree of prediction since the cont- 

20 inuous subspace between two sai^plea is now available. The pre- 
diction horison is equal to at least the size of the reference 
strip model • That means the control algorithm can be improved 
in this strip region, by performing an autoregressive control 
scheme. The operator has the visualization and command possibi- 

25 lities as shown schematically in figure 6. 

The strip is generated on the operator station using computer 
graphics tools (shadowing, material texture, translations etc) 
The strip is shown moving, in the same fashion as a video woul^ 
30 show (so called virtual re&lity) . The input to this function i 
the measured data from the measuring roll and the morphed sur- 
face • The measuring sensors touch the strip e.g- in the posi- 
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tions denoted as A/ B, C, D in figure 7. The reference model 
used for worphing has a length 1, Figure 7 shows the case when 
the prediction horizon is i. This length covers 3 measurement 
samples Ay B and The colors within the prediction area are 
mapping the difference between the actual measurement and the 
predictive values, in the usual way done with loggers , Matlab 
programs etc. For exangple red color means large error, yellow 
means medium, green is low error and the colors have gradients 
in between. The mill operator, by modifying the prediction 
horiaon in the range [0, 1] ^ and the morpbed reference model, 
can verify which case is best suited for rolling, 

To summarize, the invention involves the following? 

- a second mapping - from reference models to output space 
optimisation matrices - virtual strip is directly mapped to a 
vector space optimisation algorithm, 

- reference models enter into the computation of the main 
mapping ^ che measured flatness vector is enhanced with 
information from the virtual strip using a morphing technique, 

- use of reference pictures /models/descriptions of the strip t< 
enhance flatness error - morphing technique, 

- us© of reference pictures /models /description^ of the strip 
with an associated weight matrix to concentrate on compensatins 
zones where the problem is most significant, 

- space change. 

The method according to the invention may, at least partly, be 
performed under control of a set of computer readable 
instructions or code means contained in a computer program 
storage device for making a coi^uter or processor perform any 
of the steps of the abovedescribed method. 
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The invention may also use a computer readable product for 
carrying out the method according to the invention. 

While the present invention has been described in terms of the 
5 preferred etrfDodimentSi the invention is not limited thereto, 
but can be embodied in various ways without departing from the 
principle of the invention as defined in the claime. 
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1. Method for optimizing the measuring and determining of 
flatness faults of a strip of rolled material where a first 

5 mapping is made of the strip in at least one zone after passing 
through a mill stand, 
ebaraeterlzad Toy, 

- that a second mapping is made to optimization algorithms 
acting in actuator modal and sensor spaces using visual or 

10 image information, and 

- that the measured flatness vector is enhanced with 
information from a virtual strip, 

2. Method according to claim 1, 
15 charaeterised taF# 

that the virtual strip is mapped to a vector space optimiaatior 
algorithm. 

3. Method according to claim 1 or 2, 
20 characfcarized by, 

that a morphing techniqpie is used. 

4. Method according to any of the pjreceeding claims, 
charaeterissed by,. 

25 the result of the mapping is added 1^ morphing to the mesured 
information from a reference model. 

5. Device for optimizing the measuring and determining of 
flatness faults of a strip of rolled material after passing 

30 through a mill stand, 
eliaractarlsed hy^ 
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chat means for synchrojiizing measured data with video saitqples 
and with the currently perform optimization algoritbia is 
srranered. 

5 €. A computer program coinpriaing computer program code itieaas 
for carrying out the steps of a method according to claim 1-4. 

7. A computer readable medium comprising at least part of a 
conqputer program according co claim 6. 

0 

8. A con^ter program, according to claim 6, that is« at least 
partially, provided through a network/ such as e.g. internet. 
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^Phe invention refers to a method for optimizing the measuring 
and determining of flatness faults of a strip of rolled 
5 material where a first mapping is made of the strip in least 
one zone after passing through a mill stand • 
The invention is achieved by, 

- that a second mapping is made to optimisation algorithms 
acting in actuator model and sensor spaces using visual or 

10 image information, and 

- that the measured flatness vector is enhanced with 
information from a virtual strip. 

(Figure 3) 

15 
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